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Abstract. We present a new method for fitting peculiar velocity models 
to complete flux limited magnitude-redshifts catalogues, using the lumi- 
nosity function of the sources as a distance indicator. The method is 
characterized by its robustness. In particular, no assumptions are made 
concerning the spatial distribution of sources and their luminosity func- 
tion. Moreover the inclusion of additional observables, such for example 
the one carrying the Tully-Fisher information, is straightforward. 
As an illustration of the method, the predicted IRAS peculiar velocity 
model is herein tested using the fluxes of the IRAS 1.2 Jy sample as the 
distance indicator. The results suggest that this model, while successful 
in reproducing locally the cosmic flow, fails to describe the kinematics on 
larger scales. 



1. The method 

The application of the method is restricted to samples strictly complete up to a 
given magnitude limit mn^, i.e. the selection function in apparent magnitude is 
well described by a sharp cut-off ip{m) = 6[m\\^ — m) with 6{x) the Heaveside 
function. The probability density of the sample may be written in this case as 

dP = \h{fi, I, h) cos h dldbdfi f{M)dM 9{mii^ - m) (1) 

where ^ = m — M \s the distance modulus, h{n,l,b) the line-of-sight distribu- 
tion function, f{M) the luminosity function and A is the normalisation factor 
warranting / dP = 1. 

The milestone of the method consists in defining the random variable C, as 
follows 

c = _Zm_ , ii^-ML^M (2) 

^ F(Af|i„) ^ F(Af„„) ^ ' 

where F{M) = /j^^ f{x)dx stands for the cumulative distribution function in 
M and Miim = M\^^{^) is the maximum absolute magnitude for which a galaxy 
at distance would be visible in the sample (e.g. MiimC/u) = "'Ti-iim — A* if the 
k-correction is neglected). By definition, the random variable C, for a sampled 
galaxy belongs to the interval [0, 1]. The probability density of Eq. (||) reduces 
to 

dP = ^h{fi,l,b)F{Mnm{f^)) cosbdldbd^i x e{C)9{l-C)dC (3) 
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This equation implies the two following properties: 

• PI: (is uniformly distributed between and 1. 

• P2: C and fi are statistically independent, i.e. the distribution of C does 
not depend on the spatial position of the galaxies. 

The random variable C can be estimated without any prior knowledge of 
the cumulative luminosity function F[M). For each galaxy {Mi,^i) one can 
indeed form the subsample S'j = 5i U ^2 with Si = { (M, //) such that M < 
Mi and /U < } and ^2 = { {M,^) such that Mj < M < Mj^jj^^ and /i < /Xj }. 
By construction (see figure 1) M and /i are independent in each subsample Si. 
This implies that the following quantity is an unbiased estimate of the random 
variable C, 

C. = ^ (4) 
rii + l 

where rii is the number of objects in Sj = 5i U S2 and rj the number of objects 
in 5*1 (see Efron & Petrosian 1992). 

The radial peculiar velocity field is herein described by a linear model 
parametrized by a A^-dimensional vector (5 = {Pi,f32, ...,Pn), 

N 

Vf3 = Vf3{r) = ^Pk Vk{r) (5) 

k=l 

where vi{r), V2{r), V]y{r) is a set of functions depending on the spatial 
position r. It assumed hereafter that there exists a solution /3* fairly reproducing 
the true radial peculiar velocities of galaxies Vi = f(rj) ~ f^*(rj). For a given 
vector /?, the model dependent variables fip and Mg can be computed from the 
observed redshift z and apparent magnitude m following 

Hl3 = 5 logio ^ + 25-UI3 ; Mp = m- up ; up = -b log^g ( 1 - — ) (6) 
-Ho \ cz J 

The quantities ^p and Mp are related to the true absolute magnitude M and 
distance modulus /u via fip = n + up* — up and Mp = M — up* + up. 

It is shown in Rauzy&Hendry (hereafter RH) that any test of independence 
between (p computed following Eq. (Q) and up provides us with an unbiased 
estimate of In particular the coefficient of correlation p{Cp, up) has to vanish 
when P = f3*, i.e. 

/3 = r ^ p{Cp,up)=0 (7) 

The accuracy of this estimator can be obtained through numerical simulations by 
analysing the influence of sampling fluctuations on the coefficient of correlation. 

This estimator is clearly robust. Nothing has been assumed concerning the 
spatial distribution of sources nor on the shape of their luminosity function. In 
particular the method is free of Malmquist-like biases. Moreover the inclusion of 
a second observable parameter, for example in order to account for Tully-Fisher 
information, can be achieved without any difficulty. Finally it is worthwhile 
to mention that the method will benefit from the use of the velocity orthogo- 
nalization procedure proposed by Nusser & Davis (1995) (see RH for details). 
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Figure 1. Distance modulus versus absolute magnitude for the 60 
//m IRAS 1.2 Jy sample (5321 galaxies) 



2. Application to the IRAS 1.2 Jy sample 

The method is herein applied to the 60 ^m IRAS 1.2 Jy sample (Fisher et al. 
1995). The distance modulus versus absolute magnitude diagram is shown in 
figure 1. The peculiar velocity field model tested is the predicted IRAS velocity 
field (Strauss et al. 1992) characterized by only one parameter /3 = J^o.e/^/- 
The luminosity function of these sources does not exhibit any turnover towards 
the faint-end tail, at least within the observed range of magnitudes. Due to the 
large spread of the LF, one cannot expect very high constraints on the velocity 
model tested. However a rejection test for the P parameter can be constructed. 
We obtained that (3 > 0.7 can be rejected with a confidence level of 95% and 
P > 1.1 with a confidence level of 99%. 

In a second step, we use the observed correlation between the absolute 
magnitude M and some "colour index" defined as p = 2.5 logio(-^ioo/-^6o) (with 
FiQQ the flux at 100 iim) in order to refine the analysis. The data have been 
grouped in 8 classes by interval of p. Because of the correlation between p and M, 
the spread of the luminosity function for each of these classes taken individually 
is expected to be smaller than the spread of the global luminosity function, and 
thus the accuracy of the distance indicator improved. The random variable 
is then computed using Eq. (^ but this time class by class. The correlation 
between and the velocity modulus up is after that evaluated for the whole 
sample. 

The results are presented in figure 2 in terms of the confidence level of 
rejection for the parameter /3. The method has been first applied to the galaxies 
within the redshift range 1000-12000 km s'^ It is found that /3 G [-0.35,0.25] 
at la, and that models with (3 > 0.5 can be rejected with a confidence level of 
95%. This result is in disagreement with most of the analyses based on Tully- 
Fisher data (e.g. VELMOD on Marklll (Wilhck k Strauss 1998), ITF method 
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Figure 2. Confidence level of rejection for the parameter (3 



on SFI (Da Costa et al. 1998) favoring a value of /3 = 0.5. We interpret this 
discrepancy as follows. 

When fitting a velocity model to data, the natural weight assigned by the 
fitting procedure to each galaxy is roughly proportional to the inverse of its 
redshift. The mean effective depth of the volume where the velocity model 
is compared to data has to be estimated using these weights. For our first 
sample with z G [1000, 12000] km s~^, we find a mean effective depth 3800 km 
s~^. In order to mimic the effective volume sampled by Tully-Fisher data we 
have applied the method to a truncated version of the IRAS sample containing 
1621 galaxies with z G [500,5000] and galactic latitude \h\ > 20 (the mean 
effective depth of this sample is 2200 km s~^). Figure 2 shows that the value 
of P estimated from this truncated sample is fully consistent with the values 
obtained using Tully-Fisher data. An interpretation of these results could be 
that the predicted IRAS velocity field model, while successful in reproducing 
locally the cosmic flow, fails to describe the kinematics on larger scales. 
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